function [output] = butterworth(wp, ws, rp, rs, input)
[n, wn] = buttord(wp, ws, rp, rs);
[b, a] = butter(n, wn);
output = filter(b,a,input);
